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ABSTRACT 

We have developed a model to describe two-point correlation functions of clusters of 
galaxies in X-ray flux-limited surveys. Our model properly takes account of nonlinear 
gravitational evolution of mass fluctuations, redshift-space distortion due to linear 
peculiar velocity field and to finger-of-god, cluster abundance and bias evolution on the 
basis of the Press - Schechter theory, the light-cone effect, and the selection function 
due to the X-ray flux, temperature and luminosity limits. Applying this model in 
representative cosmological models, we have presented quantitative predictions for 
X-ray selected samples feasible from the future surveys with the X-ray satellites 
including Astro-E, Chandra, and XMM. The comparison of these predictions and the 
observed cluster clustering will place important cosmological constraints which are 
complementary to the cluster abundance and the cosmic microwave background. 

Subject headings: cosmology: theory - distance scale - dark matter - large-scale 
structure of the universe - galaxies: distances and redshifts - clusters - X-rays: galaxies 
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1. Introduction 

The abundance of cluster of galaxies is now well-established as a standard cosmological probe 
(White, Efstathiou & Frenk 1991; Jing & Fang 1994; Barbosa et al. 1996; Viana & Liddle 1996; 
Eke, Cole, & Frenk 1996; Kitayama & Suto 1997). In particular, several available X-ray catalogues 
of clusters have played an important role in placing fairly robust constraints on the cosmological 
parameters. The spatial two-point correlation function of clusters is another important target for 
cosmological researches (Bahcall & Soneira 1983; Klypin & Kopylov 1983; Bahcall 1988; Bahcall & 
Cen 1993; Ueda, Itoh & Suto 1993; Jing et al. 1993; Watanabe, Matsubara & Suto 1994; Borgani, 
Plionis, & Kolokotronis 1999). In fact the idea of biased galaxy formation by Kaiser (1984) was 
proposed originally to reconcile the stronger spatial correlation of clusters with those of galaxies. 
The previous studies in this context, however, have been mainly based on the optically selected 
samples, which are likely to be contaminated by the projection effect and thus the selection 
function of which is difficult to evaluate precisely. 

The X-ray flux-limited catalogues of clusters are now becoming available with ROSAT (e.g. 
Ebeling et al. 1997; Rosati et al. 1998) and are expected to increase in their sample volume in 
near future with the X-ray satellites including Astro-E, Chandra, and XMM. Such well-controlled 
catalogues are ideal to revisit the cluster correlation functions with unprecedented precision. 
The proper comparison with such data, however, requires better theoretical predictions which 
take account of the selection function of X-ray clusters (Kitayama, Sasaki & Suto 1998), the 
luminosity and time dependent bias (Mo & White 1996; Jing 1998; Moscardini et al. 1998), 
the light-cone effect (Matarrese et al. 1997; Matsubara, Suto &; Szapudi 1997; Nakamura, 
Matsubara & Suto 1998; Yamamoto & Suto 1999; Moscardini et al. 1999a) and the redshift-space 
distortion (Hamilton 1998; Matsubara & Suto 1996; Suto et al. 1999; Nishioka & Yamamoto 1999; 
Yamamoto, Nishioka & Suto 1999; Magira, Jing & Suto 2000). The theoretical prescriptions for 
those effects have been developed and become available recently. This motivates us to present 
detailed predictions for the two-point correlation functions for X-ray flux-limited samples of 
clusters of galaxies combining the cluster abundances in a fully consistent fashion. 

Moscardini et al. (1999a) recently performed a similar work, and our present paper differs 
from theirs in several aspects. First, we take account of the linear and nonlinear redshift-space 
distortion due to the peculiar velocity field. Second we adopt a formula for the light-cone effect 
derived by Yamamoto & Suto (1999) which is the one-dimensional integration over the redshift 
in the cluster sample. Finally, we are mainly interested in future surveys which probe a higher 
redshift (z ~ 1) universe, and therefore focus on a range of the X-ray flux limit a few magnitudes 
fainter than they considered. 

Koyama, Soda & Taruya (1999) suggested a presence of primordial non-Gaussianity from 
the combined analysis of the cosmic microwave background, the abundance of X-ray clusters at 
z = and 0.3, and the correlation length of optical cluster samples. While this interpretation is 
interesting, the definite conclusion requires the proper account of the selection effect as well as the 
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theoretical modeling and the statistical limitation of the available sample. In this sense, optically 
selected cluster samples currently available are still far from satisfactory, and the analysis of the 
upcoming X-ray selected catalogues is essential. 



2. Modelling the two-point correlation functions of X-ray clusters 

2.1. Linear and nonlinear redshift-space distortion 

The observable two-point correlation functions are inevitably distorted due to the presence of 
the peculiar velocity field. We take into account this redshift-space distortion following Cole et al. 
(1994), Magira et al.(2000) and Yamamoto et al. (1999). 

Our key assumption is that the bias of the cluster density field relative to the mass density 
field is linear and scale-independent: 

5 c i(x, z) = b c \(z) <5 mass (x, z). (1) 

In this case, the power spectrum of the corresponding cluster samples in redshift space is well 
approximated as 

l 2 



1 + {k^a v ) 2 /2 



[bd(z)] 2 I*La(k,z), (2) 



where \i the direction cosine of the wavenumber vector and the line-of-sight of the fiducial observer, 
and -Pmass(^' z ) i s the mass power spectrum in real space. The numerator in equation @/ expresses 
the linear redshift-space distortion (Kaiser 1987), where (5 is defined by 

b c i{z) (ilna(z) 

and D\[z) is the linear growth factor normalized to be unity at present. The denominator in 
equation @ takes account of the nonlinear redshift-space distortion {finger- of- God) assuming 
that the one-point distribution function of the peculiar velocity is exponential with the velocity 
dispersion of a v (z). 

Averaging equation (pi) over the angle with respect to the line-of-sight of the observer yields 
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with k(z) = ka v {z)/ Hq. Finally the corresponding two-point correlation function of clusters in 
redshift space is computed as 

1 f°° 

^(R, z ) = —J^ dkk 2 P! l (k,z)j (kR), (8) 

where jo(kR) is the spherical Bessel function. In what follows, we adopt nonlinear evolution of 
the mass fluctuations using the Peacock & Dodds (1996; PD) fitting formula for P^ ass (/i;, z) unless 
otherwise stated. 



2.2. Temperature and luminosity of X-ray clusters 

The X-ray luminosity function of galaxy clusters at z ~ is determined to a good precision 
from existing observational catalogues (Burns et al. 1996; Ebeling et al. 1997). In fact, the 
Press-Schechter theory applied to the CDM models reproduces quite well the X-ray luminosity 
function, as well as the X-ray temperature function, provided that the amplitude of the mass 
fluctuation at 8/i _1 Mpc, <7g, is related to the density parameter Qq and the cosmological constant 
Ao as follows (Kitayama & Suto 1997; see also Viana & Liddle 1996; Eke, Cole & Frenk 1996): 

o -0.35-0.82n +0.55nj? 1 „ A 

fJ 8 -U.54x<; 0.28-0.91^+0.68^ , , ■ W 

U Q (A = UJ 

In this paper, we follow the prescription of Kitayama & Suto (1997) to relate the mass M of 
the "Press-Schechter objects" to the X-ray luminosity L (bolometric or band-limited) of the real 
clusters of galaxies. Although the one-to-one correspondence of those two species is a non-trivial 
assumption, it is regarded as a fairly successful approximation, at least empirically. We first relate 
the total mass M of the dark halo of a cluster to the temperature of hot gas, T gas , assuming the 
virial equilibrium: 

lim p GM 



3r vir (M, z)' 



= 5-27(1 +z) ( ^) 1/3 (-M—) ^ ^ 3 keV. (10) 



Vl8vrV VIO 15 ^- 1 ^©. 

where k B is the Boltzmann constant, G is the gravitational constant, m p is the proton mass, fi is 
the mean molecular weight (we adopt fi = 0.59), and 7 is a fudge factor of order unity (we adopt 
7 = 1.2). In the above, z is the redshift of the cluster which we assume is equal to the cluster 
formation epoch for simplicity (see Kitayama & Suto 1996 for more discussion on this point). 
While the cluster temperature may not be isothermal, it is not important in our analysis. In fact, 
the above temperature should be interpreted as an average temperature of the cluster. In fact, 
our latest hydrodynamical simulations (Yoshikawa, Jing & Suto 1999) show that the mass- and 
emission-weighted temperatures of simulated clusters satisfy the above relation with 7 = 1.2 and 
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1.6, respectively. Further details of the temperature-mass relation and its non-isothermal effect 
are discussed in Makino, Sasaki & Suto (1998), Yoshikawa, Itoh, & Suto (1998), Suto, Sasaki & 
Makino (1999) and Yoshikawa & Suto (1999). We compute A v ; r , the ratio of the mean cluster 
density to the mean density of the universe at that epoch using the formulae for the spherical 
collapse model presented in Kitayama & Suto (1996). 

Next we transform the temperature to the luminosity of clusters using the observed luminosity 
- temperature relation: 

L bol = L 44 (^) Q (l + ^ lO 44 ^" 2 erg sec- 1 . (11) 

While a simple theory on the basis of the self-similar cluster evolution predicts the slope a = 2, we 
adopt L44 = 2.9, a = 3.4 and £ = on the basis of recent observational indications (e.g., David et 
al. 1993; Ebeling et al. 1996; Ponman et al. 1996; Mushotzky & Scharf 1997). Then we translate 
L hol (T g , s ) into the band-limited luminosity L band [T gas , E u E 2 ] as 

£ b a nd [T gas , E a (l + z), E b {l + z)} = L bol (T gas ) x /[T gas , E a (l + z), E b (l + z)}, (12) 

where f[T gas , E±, E 2 ] is the band correction factor which takes account of metal line emissions 
(Masai 1984) in addition to the thermal bremsstrahlung. 

Finally the source luminosity -L band at z is converted to the observed flux S in an X-ray energy 
band [E a ,E b \. 

S[Ea > Eb] = A^aW) (13) 

where d L (z) is the luminosity distance. Throughout the present paper, we use the 0.5 — 2.0keV 
band for the flux assuming the abundance of intracluster gas as 0.3 times the solar value. The 
luminosity distance is related to the comoving distance d c and the angular diameter distance d A 
as d L {z) = (1 + z)d c {z) = (1 + z) 2 d A (z); 



d c (z) = < 
where r{z) is the radial distance: 

and H (z) the Hubble parameter at redshift z: 



' Hq 1 sin {H r^n + A - l)/V^o + A - 1 (fio + A > 1) 

r (ft + A = l) , (14) 

k Hq 1 sinh {H ry/1 - Sl - Aq)/\/1 - ^0 - A (Oq + A < 1) 



r(z) = '„ wr <15> 



H(z) = H Jn (l + zf + (1 - n - A )(l + zf + A . (16) 



In what follows, we mainly consider three representative models; SCDM (Standard CDM) 
with (n ,X ,h,a 8 ) = (1.0,0.0,0.5,0.56), LCDM (Lambda CDM) with (0.3,0.7,0.7,1.04), and 
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OCDM (Open CDM) with (0.45,0.0,0.7,0.83), where h is the Hubble constant H in units of 
lOOkm/s/Mpc. It should be noted that equation (^) is valid for h = 0.7, strictly speaking. For 
that reason, we adopt = 0.56, instead of 0.54, to match the cluster abundance for SCDM with 
h = 0.5. 

The relations among mass, temperature, luminosity and flux of X-ray clusters which we 
describe in the above do not depend on the power spectrum. For defmiteness, however, we choose 
the same cosmological parameters as the above three models, and plot the results. Figure [l 
shows the temperature Tx and the bolometric luminosity Lt>oi> as a function of cluster mass M 
at z = 0.01 (dashed lines), 0.5 (dotted lines) and 1.0 (solid lines). Figure ^ plots Tx and M as 
a function of X-ray flux So.5-2.0 ^ or dusters at z = 0.01 (dashed lines), 0.5 (dotted lines) and 
1.0 (solid lines). Figure || shows Tx, Tbol and M of clusters at z with the observed X-ray flux 
So.5-2.0 = 10~ 13 , 10~ 14 and 10 _15 erg/s/cm 2 in solid, dotted and dashed curves, respectively. 



2.3. Evolution of luminosity dependent bias 



In order to complete a model for £fi(i2, z; > S\- im ), the two-point correlation function of X-ray 
clusters on a constant-time hypersurface at z, one has to specify a model for bias. Here we adopt a 
fitting formula for bias of Jing (1998), which improves an original proposal by Mo & White (1996) 
on the basis of high-resolution N-body simulations: 

0.06-0.02n eff 



b(z,M) 



0.55* 



A 4 (M,z 



+ 1 



1 



1 



+ 



A 2 (M, 



(17) 



In the above, A(M, z) is the density fluctuation smoothed over a (top-hat) mass scale of M at the 
redshift z, n e ff is the effective slope of the linear density power spectrum at the mass scale M 
(see Jing 1998), 5 C = 1-69 is the critical density contrast for the spherical collapse. To be specific, 
Jing (1998) carried out several simulations in four scale- free models and three representative CDM 
models employing 256 3 particles. Each model is simulated with three or four different realizations. 
The dark matter halos are identified using the friends-of-friends algorithm using the linking length 
of 0.2 times the mean particle separation. Thus the above formula ( |l7|) applies for the virialized 
dark halo defined according to the Press - Schechter theory, and is appropriate for the clusters of 
galaxies. 

Combining the results in the previous subsection, we translate the above bias factor into a 
function of X-ray flux limit according to 



b cS (z, > Sum) 
n PS (z,M) 



M lim (5 lim ) 



dM b(z,M) n PS (z,M) 



dM n PS {z,M) 



2X 1 / 2 3n H% 5 C dA(M,z) 
8itGM A 2 {M,z) dM 



cxp 



2A 2 (M, z) 



(18) 
(19) 
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where n PS (z, M) is the Press - Schechter mass function. Given the flux limit S\i m , the corresponding 
luminosity L^ m (z) for a cluster located at z is computed from equation ( |l3| ) with the luminosity 
distance di,(z). Then we solve equations (|TT|) and (12) for Ti im (z) = T gas , and finally obtain 



M\\ m (z) using equation (10) 



Figure || plots the evolution of b e g(z, > M), b(z,M), b e g(z, > S\i m ) and b(z, Sn m ) for our 
model of X-ray clusters. As noted in Figure 3 of Jing (1998), the fitting formula is accurate within 
10 percent for M £ 10 14 h~ l M & on linear scales. We checked the accuracy of the fit using the 
more recent numerical simulations, and made sure that the fit is accurate better than 20 percent 
even at M = 1O 15 /i _1 M . In any case, the fraction of such massive clusters is significantly smaller 
according to equation fll9| ) , and the difference of the fit and simulations to that level hardly changes 
our results in reality. We should note here, however, that the bias factor becomes scale-dependent 
for scales below ~ 5/i _1 Mpc, and then our assumption of linear and scale-independent bias breaks 
down. Thus while our predictions on linear scales ^ 5/i _1 Mpc, which we are mostly interested 
in, are reliable, those on smaller scales should be interpreted simply to show the various effects 
qualitatively. 



2.4. The light-cone effect 

The two-point correlation function on the light-cone is properly formulated by Yamamoto & 
Suto (1999). In the present context, their formula yields the following expression for the two-point 
correlation functions of clusters brighter than the X-ray flux- limit Sn m : 

~d~Vc 



&°- c i(R; > s hm ) = ^ — d " (20) 



r^uKz) 



where R is the comoving separation of a pair of clusters, z max and z m i n denote the redshift 
range of the survey, and (R, 

z ) > "Slim) is the corresponding two-point correlation function on 
a constant-time hypersurface at z in redshift space (eq.Q). The comoving number density of 
clusters in the flux-limited survey, no(z; > Sn m ), is computed by integrating the Press - Schechter 
mass function as 



n (z;> S hm )= n PS (M,z)dM. (21) 

JM nm (S lim ) 



Finally the comoving volume element per unit solid angle is 

dV c _ d 2 c (z) 
~dz ~ H{z) 



dv c _ dUz) (22) 



While equation ( |20| ) was originally derived for the cosmological models with flat spatial 
geometry (Yamamoto & Suto 1999), we assume that the formula is valid for the cosmological 
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models with open geometry. Note also that this expression for the light-cone effect looks rather 
different from that adopted by Mataresse et al. (1997) and Moscardini et al. (1999a), but 
that both lead to fairly similar results quantitatively; see Yamamoto & Suto (1999) for detailed 
comparison and discussion. 

The formula (20) is useful in making theoretical predictions, but no(z) itself is not directly 
observable (unless the cosmological parameters are specified). Instead it can be rewritten in terms 
of the redshift distribution of clusters per unit solid angle with an X-ray flux S > S\i m : 

f^ H(z)dz fdN\\ s 

thC (Tf.^O \ _ Jz mi n d 2 c {z) \dz J , , 

SX-cll-Ki > ^limj — — —, : . . — o ) \ 16 ) 



r z ™* H(z)dz (dN\ 2 



and dN/dz which is related to uq(z) as 



dN dV c . A , 

-j^-{z,> Sn m ) = n (z; > 6 Um ) — . (24) 

Figure || plots the redshift distribution function for X-ray flux-limited clusters with 
<Siim = 10~ 13 , 10~ 14 and 10 -15 erg/s/cm 2 . It is interesting to note that My im (z) corresponding to 
a fixed Sum decreases for z ^ 1. In fact this behavior is understood from equations (10) to (|i~3|); 

M lim (z) oc T^(z)/(1 + zfl* oc L?/ m 2Q (z)/(l + z )*l* oc d 3 L /a (z)/(l + zfl\ 



3. Predictions for correlations in X-ray flux-limited surveys of clusters of galaxies 

3.1. Nonlinear evolution, redshift-space distortion and light-cone effect 

While our predictions based on the modeling described in the previous section include various 
important effects, it would be instructive to discuss them separately first. For that purpose, 
we plot in Figure ^ several different predictions for two-point correlation functions in which 
some of the effects are artificially turned off; linear and nonlinear mass correlations in real space 
at z = using the Bardeen et al. (1986; BBKS) and PD formulae for mass power spectra, 
cluster correlations with linear redshift-space distortion (Kaiser 1987) and with full redshift-space 
distortion, at z = using the fitting formula of Mo, Jing & Borner (1997) to compute a v . These 
should be compared with our final predictions on the light-cone in redshift space (eq.||). 

The qualitative features illustrated in Figure ^ can be understood as the combination of 
the following effects^]; the nonlinear gravitational evolution increases the correlation on small 
scales, while redshift-space distortion decreases (enhances) the amplitudes on small (large) scales. 



As mentioned earlier, the behavior of the correlation functions on nonlinear scales should not be trusted because 
we have not properly taken into account the possible scale-dependence of the bias below the scales. 
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The correlation length defined by eq.[^| below, for inscance, is enhanced in redshift space by 
50 percent for SCDM and 20 percent for LCDM and OCDM in comparison with those in real 
space. The light-cone effect averages the amplitude over a range of redshift which is generally 
expected to decrease the correlation if the clustering amplitude at higher z decreases according to 
linear theory. In reality, however, b e s(z) increases more rapidly than the linear growth rate Dy(z) 
at higher z, and therefore the clustering amplitude evaluated on the light cone for a given Sii m 
increases as z (see Fig.|| below). 

Nevertheless this is not what we observe in a flux-limited survey. If the survey flux-limit Sii m 
becomes smaller, the sample includes both less luminous clusters at low z and more luminous 
clusters at high z. While the latter shows stronger bias, the former should exhibit weaker bias. 
Thus averaging over the light-cone volume, their effect on the overall clustering amplitude is fairly 
compensated. This explains why the overall results are not so sensitive to S\i m as shown in left 
panels of Figure [7|. 



3.2. Effect of the selection function 

If the temperature of an individual cluster in an X-ray flux-limited sample is determined 
observationally, one can construct a temperature-limited subsample. Similarly one can construct 
a luminosity-limited subsample from the flux-limited sample with the redshift information of 
each cluster. If the underlying luminosity-temperature relation (eq. [|llj) had no dispersion, 
these two subsamples would be essentially the same. In reality, however, the finite amount of the 
dispersions would lead to different predictions for those subsamples, which will provide further 
and independent information on the cluster models. This is particularly the case here since the 
bias is highly sensitive to the mass, and therefore to the temperature and luminosity, of a cluster. 

Figure |?] shows our predictions for £ L C ^ (R) for cluster samples selected with different flux-limit 
"Slim {left panels), and with temperature and absolute bolometric luminosity limits, T\\ m and 
L\\m {middle and right panels). For the latter two cases, Sy im = 10~ 14 erg/s/cm 2 is assumed for 
definiteness. As explained in the previous subsection, the results are insensitive to 5i; m , but very 
sensitive to T\ im and L\ im , reflecting the strong dependence of the bias on the latter quantities. 

To see this in a somewhat different manner, we plot in Figure || the correlation length r c o(Sn m ) 
defined through 

^ C (r c0 ;>5 lim ) = 1, (25) 

as a function of Sn m , T\ im and L\ im . For the latter two cases, Sn m = 10 _14 erg/s/cm 2 is assumed as 
in Figure [7|. Again it is clear that the results are fairly insensitive to Sy im , but are sensitive to the 
bias factor, which results in the strong dependence of r C Q on T\\ m and L\\ m . 

The value of r c o(Sn m ) also depends on the depth of the survey via the light-cone effect and 
the evolution of bias and mass fluctuations (Fig.||). As a result of the several competing effects, 
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fco(Siim) increases, albeit very weakly, as z max becomes larger. 



3.3. Dependence on Sl 

For a cosmological application of the present result, it is interesting to examine how the 



r c,o(5'iim) depends on S7o- This is summarized in Figure 10, where we fix the value of the 
fluctuation amplitude as adopting the cluster abundance constraint (eq.|J), and consider both 
Aq = 1 — SIq (thick lines) and Ao = (thin lines). We set the shape parameter of the spectrum 
T as Qohexp[— f?b(l + V^h^-Q 1 )] with fl^h 2 = 0.015 and h = 0.7. If we adopt such spectra, the 
correlation length in Ao = 1 — model is generally larger than that in Ao = model, largely 
reflecting the dependence of ug on Ao (eq. ||). Again the results are not sensitive to the flux limit 
Siim. The dependence on Qq is rather strong, and these predictions combined with the future 
observational results will be able to break the degeneracy of the cosmological parameters. 



4. Conclusions and discussion 

We have presented a detailed methodology to predict the two-point correlation functions 
for X-ray flux-limited samples of clusters of galaxies, fully taking into account the redshift-space 
distortion, nonlinear gravitational evolution of mass fluctuations, evolution of bias, the light-cone 
effect and the observational selection function. While our method is similar to the recent model of 
Moscardini et al. (1999a) in many respects, the most important difference is that we incorporated 
both linear and nonlinear redshift-space distortion following the prescription of Suto et al. (1999), 
Yamamoto, Nishioka & Suto (1999) and Magira, Jing & Suto (2000). 

The predictions for correlations of clusters are most sensitive to the model of bias among 
others. Fortunately, the Press - Schechter theory provides a quite reliable bias model for clusters 
(Mo & White 1996; Jing 1998), which is in marked contrast with a model for galaxy and 
quasar bias. Thus one can make quantitative and detailed model predictions once a set of the 
cosmological parameters are specified. These predictions can be checked against the flux-limited 
samples of clusters from the future X-ray satellites including Astro-E, Chandra, and XMM. It is 
particularly interesting to probe the value of f^o and the degree of non-Gaussianity by combining 
this comparison and the cluster abundance (Kitayama, Sasaki, & Suto 1998; Koyama, Soda, & 
Taruya 1999). 

We should note, however, that the current methodology is reliable only for linear scales; on 
small scales ( £ 5/t Mpc), the approximation of the linear and scale- independent bias breaks 
down. In addition, the projection effect becomes inevitably important below a few h~ 1 Mpc scales 
even for the X-ray selected samples, and the observational data analysis itself becomes non-trivial. 

While the selection function for the future survey would not exactly match those in our 
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examples, it is straightforward to make suitable predictions on the basis of the present formalism 
at least on linear scales, and we hope that we have already presented all the basic and qualitative 
features of the clustering statistics. 

After submitting this paper, Moscardini et al. (1999b) posted a quite similar paper on 
the clustering of X-ray selected galaxy clusters. They predicted the correlation lengths at 
Sum = 10~ 14 erg/s/cm 2 of 7 and 12 /i _1 Mpc in SCDM and LCDM models, respectively (their 
Fig. 7), which should be compared with our predictions of 9 and 17 h~ 1 Mpc (Fig.|||). While 
our models adopt slightly different values for the shape parameter T, and the normalization 
of fluctuations, erg, it is clear that our correlation lengths are systematically larger by 30 ~ 40 
percent level. This should be ascribed to be the redshift-distortion effect which they neglect in 
their analysis. In fact, this difference is comparable to the cosmological model dependence which 
they showed in their Figure 7, clearly suggesting the importance of the effect to this level. Apart 
from this difference, however, our results are consistent with each other. 

We thank an anonymous referee for detailed comments which improved the presentation of the 
earlier manuscript. Y.P.J, and T.K. gratefully acknowledge the fellowship from the Japan Society 
for the Promotion of Science. Part of numerical computations was carried out on VPP300/16R 
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Observatory, Japan, as well as at RESCEU (Research Center for the Early Universe, University 
of Tokyo) and KEK (High Energy Accelerator Research Organization, Japan). This research was 
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Fig. 1. — The X-ray temperature and luminosity as a function of cluster mass for three 
representative sets of cosmological parameters; (f2o,Ao) = (1.0,0.0), (0.3,0.7), and (0.45,0.0) from 
left to right panels. The X-ray temperature (Upper panels) and the bolometric X-ray luminosity 
(Lower panels) are plotted at z = 1.0 (solid lines), 0.5 (dotted lines) and 0.01 (dashed lines). 
We adopt the ratio of gas to virial temperature 7 = 1.2 (eq. fTo|), and the non-evolving (( = 0) 
luminosity-temperature relation (eq.|ll|) with the amplitude L44 = 2.9 and the power-law index 
a = 3.4. 
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Fig. 2. — The X-ray temperature and cluster mass as a function of X-ray flux in the 0.5 — 2 keV 
band. Different lines indicate results at different redshifts; z = 1.0 (solid), 0.5 (dotted) and 0.01 
(dashed). 
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Fig. 3. — Temperature, bolometric luminosity and mass of clusters corresponding to the X-ray flux- 
limit Sii m as a function of redshift z; Sn m = 1CT 13 (solid lines), 10~ 14 (dotted) and 10 _15 erg/s/cm 2 
(dashed). 
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Fig. 4. — Evolution of bias for clusters selected by mass and X-ray flux in SCDM, LCDM and 
OCDM models; Upper panels: bias for different mass limits, Mam = 10 14 fc -1 Mo, lO 13 fc- 1 Af 
and 10 12 /i _1 Mq in solid, dotted and dashed lines, respectively. Thin lines indicate the result for 
b(z, Miim) while thick lines for the effective bias b e s(z,> Mu m ). Lower panels: bias for different 
flux limits, Sii m = 10~ 13 , 10 -14 and 10 _15 erg/s/cm 2 in solid, dotted and dashed lines, respectively. 
Thin lines indicate the result for b(z, Su m ) while thick lines for the effective bias b e s(z, > Sy m ). 
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Fig. 5. — Redshift distribution functions for X-ray flux-limited clusters; Sn m = 10 , 10 and 
10~ 15 erg/s/cm 2 in solid, dotted and dashed lines, respectively. 




Fig. 6. — Light-cone and redshift-space distortion effects on two-point correlation functions 
of clusters; real-space mass correlation functions at z = in linear theory using the BBKS 
transfer function (dashed lines) and in a nonlinear model by PD (dashed lines). Redshift-space 
mass correlation functions at z = with the Kaiser distortion (dot-dashed lines) and with the 
Kaiser distortion and finger-of-god (thin solid lines). Thick solid lines indicate our predictions 
for the cluster correlation functions on the light cone in redshift space (the X-ray flux limit 
Snm = 10- 14 erg/s/cm 2 ). 
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Fig. 7. — Two-point correlation functions of clusters in SCDM (Top panels), LCDM (Middle 
panels), and OCDM (Bottom panels) for different selection criteria. Left panels: the X-ray flux 
limit Si im = 1CF 13 (solid lines), 10~ 14 (dotted) and 10 _15 erg/s/cm 2 (dashed). Central panels: 
clusters with the temperature larger than 1 (solid), 3 (dotted) and 6keV (dashed) in the X-ray flux- 
limited sample (S\\ m = 10 _14 erg/s/cm 2 ). Right panels: clusters with the bolometric luminosity 
larger than 10 45 (solid), 10 (dotted) and 10 43 /i _2 erg/s/cm 2 (dashed) in the X-ray flux- limited 
sample (Sy im = 10~ 14 erg/s/cm 2 ). 
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Fig. 8. — Correlation lengths of clusters r c o in SCDM (Top panels), LCDM (Middle panels), and 
OCDM (Bottom panels) for different selection criteria. Left panels: r C Q for the X-ray flux-limited 
clusters. Dashed lines refer to the results in which the bias parameter is set to unity, while solid 
lines indicate our predictions. Central panels: r C Q as a function of the temperature limit for the X- 
ray flux-limited sample of Sy im = 10~ 13 (solid lines), 10 -14 (dotted) and 10 _15 erg/s/cm 2 (dashed). 
Right panels: r C Q as a function of the bolometric luminosity limit for the X-ray flux-limited sample 
of S hm = 1(T 13 (solid lines), 10" 14 (dotted) and l(T 15 erg/s/cm 2 (dashed). 
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Fig. 9. — Correlation lengths of clusters in SCDM (Left panels), LCDM (Middle panels), and 
OCDM (Right panels) as a function of the survey depth, z max . The X-ray flux-limit S\- im is 10~ 13 
(solid), 10~ 14 (dotted) and 10 _15 erg/s/cm 2 (dashed lines). Thin lines indicate the predictions in 
the b(z) = 1 model for comparison. 




Fig. 10. — Correlation lengths of clusters as a function of f^o- The shape parameter of the spectrum 
r is fixed as £lohexp[— £1^(1 + \/2M7q 1 )] with Q,\ 3 h? = 0.015 and h = 0.7. The X-ray flux-limit S\ ira 
is 10 -13 (solid lines), 10~ 14 (dotted) and 10 _15 erg/s/cm 2 (dashed). For each Sy im , we plot the case 
of Aq = 1 — hi thick lines, and Aq = in thin lines. 



